SCALAR, VECTORIAL AND TENSORIAL DAMAGE PARAMETERS FROM 
THE MESOSCOPIC BACKGROUND 

C. PAPENFUSSi, P. VAN2 

Abstract. In the mesoscopic theory a distribution of different crack sizes and crack 
orientations is introduced. A scalar damage parameter, a second order damage tensor and 
a vectorial damage parameter are defined in terms of this distribution function. As an 
example of a constitutive quantity the free energy density is given as a function of the 
damage tensor This equation is reduced in the uniaxial case to a function of the damage 
vector and in case of a special geometry to a function of the scalar damage parameter. 



1. Introduction 

1.1. Phenomenological definitions of damage parameters. Numerous damage models 
have incorporated scalar, vectorial or tensorial damage variables that can be characterized 
at the macro-scale, for example, by the change in compliance. A scalar damage parameter 
has been introduced lU |2] [3] |4l to account for the decrease in the stiffness of the ma- 
terial with progressing damage. Two scalar damage parameters have been proposed f5\ 
to account independently for the change in hydrostatic energy and the remaining part of 
the elastic energy with increasing damage. A different reason for introducing two scalar 
damage parameters was in \6i to account for healing of cracks under compression. In com- 
posites and fiber reinforced materials it is reasonable to introduce independent scalar order 
parameters for the prescribed directions, given by the fiber orientation. 

A second order damage tensor has been defined Q, accounting for the reduction of the 
effective surface area, which transmits forces. The resulting effective stress is expressed 
in terms of the damage tensor |7|. For definitions of a second order damage tensor see 
also m m m \W\ [TTI . For parallel microcracks a second order damage tensor has been 
associated with the dyadic product of crack orientation with itself times a scalar parameter 
lfT2l . This definition coincides with our definition from the mesoscopic point of view in 
the special case of parallel microcracks. 
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A fourth order damage tensor has been introduced. It can be understood as mapping 
the elastic tensor of the virgin material to the elastic tensor of the damaged material, or as 
mapping the respective stress tensors. For a summary of damage parameters of different 
orders see also ifTSl fT4ll . 

For a constitutive theory of damaged materials with a thermodynamic background see 
lfT5l[T6l . A thermodynamic theory of damage, including the interpretation of failure as loss 
of thermodynamic stability, can be found in ifTTll . For a comparison to experimental results 
see ifTSl . 

An alternative choice of damage variable is one that incorporates salient aspects of dam- 
age morphology in its definition. Such "micro-mechanically- inspired" damage models in- 
volving scalar, tensor, or "Fabric tensor" representations of damage have been introduced 
in the study of heterogeneous materials containing voids or various crack-Uke surface dis- 
continuities lfT9ll20ll2ni22ll23ll24l. 

Our aim here is to show how damage parameters of different tensor order can be defined 
from the mesoscopic background. The different damage parameters correspond to different 
levels of macroscopic approximation of the mesoscopic distribution of crack sizes and 
orientations. As it has been shown in Il25ll26l in the case of a scalar damage parameter, 
the mesoscopic theory leads not only to the definition of damage parameters, but also to 
equations of motion for them. On the example of the free energy density we will show the 
general form of a constitutive equation for the different choices of a damage parameter. In 
the case of a rotation symmetry of crack orientations, the different forms of constitutive 
equation can be reduced to a form with two scalar parameters: the average crack size and 
a scalar orientational order parameter. 

1.2. Mesoscopic theory of complex materials and application to material damage. 

The mesoscopic theory has been developed in order to deal with complex materials within 
continuum mechanics ll27l . The idea is to enlarge the domain of the field quantities by 
an additional variable, characterizing the internal degree of freedom connected with the 
internal structure of the material. In a simple model the micro-crack is described as a flat, 
rotation symmetric surface, a so called penny shaped crack. In addition we make here the 
following simplifying assumptions: 

(1) The diameter of the cracks is much smaller than the linear dimension of the con- 
tinuum element. Under this assumption the cracks can be treated as an internal 
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Structure of the continuum element. The cracks are assumed small enough that 
there is a whole distribution of crack sizes and orientations in the volume element. 

(2) The cracks are fixed to the material. Therefore their motion is coupled to the 
motion of representative volume elements. 

(3) The cracks cannot rotate independently of the material, i. e. the rotation velocity 
is determined by the antisymmetric part of the time derivative of the deformation 
gradient of the surrounding material. 

(4) The number of cracks is fixed, there is no production of cracks, but very short 
cracks are preexisting in the virgin material. 

(5) The cracks cannot decrease area, but can only enlarge, meaning that cracks cannot 
heal. 

To summarize our model the micro-crack is characterized by a unit vector n representing 
the orientation of the surface normal and by the radius I of the spherical crack surface. 
These parameters will be taken as the additional variables in the mesoscopic theory. 

Beyond the use of additional variables the mesoscopic concept introduces a statistical 
element, the so-called mesoscopic distribution function. In our case this is a distribution 
of crack lengths and orientations in the continuum element at position x and time t, called 
here crack distribution function (CDF). The distribution function is the probability density 
of finding a crack of length I and orientation n in the continuum element. The elements are 
material elements, including the same material and the same cracks for all times. Macro- 
scopic quantities are calculated from mesoscopic ones as averages over crack sizes and 
crack orientations. 

1.3. Mesoscopic balance equations. Field quantities such as mass density, momentum 
density, angular momentum density, and energy density are defined on the mesoscopic 
space. For distinguishing these fields from the macroscopic ones we add the word "meso- 
scopic". In addition to mass density we introduce the crack number density N as the 
density of an extensive quantity. The mesoscopic crack number density N{1, n, x, t) is the 
number density, counting only cracks of length I and orientation n. 

Balance of crack number 
In our model the cracks move together with the material element. Therefore their flux is 
the convective flux, having a part in position space, a part in orientation space, and a part 
in the length interval. There is no production and no supply of crack number. Therefore 
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we have for the crack number density N: 

(1.1) liV(.) + V..{7V(.)^(a.,t)} + V„.{iV(.)«(-)} + ^|(^'/>(-)) = 0. 

We have used spherical coordinates for the mesoscopic variables crack length I £ [0, oo] 
and crack orientation n € S'^, and we represent the divergence with respect to the meso- 
scopic variables in spherical coordinates. V„ denotes the covariant derivative on the unit 
sphere, v is the material velocity. In our model all cracks within the continuum element 
move with this velocity, it(-) = n is the orientation change velocity, which is not the same 
for all cracks in the continuum element. It is related to the angular velocity Lj{x,t) by the 
relation 

(1.2) m(-) = u: X n 

This angular velocity is the same for all cracks in the element. It is determined by the 
rotation of the surrounding material. 

1 .4. Definition of the distribution function and equation of motion. Due to its defini- 
tion as probability density the distribution function is the number fraction 

n ■2\ fn +^ N{l,n,x,t) 

(1.3) /(/,n,.,t) = ^^^-^ , 

in volume elements, where the number density N{x, t) is non-zero. Here N{x, t) is the 
macroscopic number density of cracks of any length and orientation. As the distribution 
function in equation (11.3b is not well defined if N{x, t) = 0, we define in addition that in 
this case /(/, n, x, t) — 0. As there is no creation of cracks in our model the distribution 
function will be zero for all times in these volume elements. In all other volume elements 
with a nonzero crack number it is normalized 

(1.4) / / f{l,n,x,t)l^cfndl^l . 
Jo Js^ 

With respect to crack length it is supposed that the distribution function has a compact 
support, meaning that in a sample there cannot exist cracks larger than the sample size. 

We obtain from the mesoscopic balance of crack number density a balance of the CDF 
f{l, n, x, t), by inserting its definition: 
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— /(/, n, X, t) + V, ■ {v{x, t)f{l, n, x, t)) + 
ot 

(1.5) V„ • iu{x,t)f{l,n,x,t)) + {lHf{l,n,x,t)) = 0. 

The right hand side is equal to zero, as for the co-moving observer the total number of 
cracks in a volume element does not change in time. 

A growth law for the single crack I is needed in equation ( 11.51 ). For example the Rice- 
Griffith dynamics, which is motivated from macroscopic thermodynamic considerations, 
has been applied. 

In ||28l the mesoscopic theory has been specialized to damaged material with penny 
shaped cracks. The balance equations and the differential equation for the crack size distri- 
bution function have been derived. With Rice-Griffith differential equation for the size of 
a single crack, the time evolution of the whole distribution of cracks under load has been 
investigated, as well as the evolution of the average crack size ll29l . In ll30l two different 
growth laws for the single crack under load have been considered. Finally, the dynamics of 
a second order damage tensor has been derived in ||26l under the assumption of a simplified 
single crack growth law under an effective stress. 

2. Mesoscopic definitions of damage parameters of different orders 

Definitions of a scalar damage parameter, a vectorial parameter and a damage tensor 
are given, based on the mesoscopic distribution function. A scalar damage parameter is 
the average crack length. The crack growth introduces an anisotropy into the material. In 
order to describe this anisotropy, it is necessary to define a vectorial or a tensorial damage 
parameter. Starting from the mesoscopic distribution function, the more natural way is to 
define a damage tensor of second order. The tensor character is introduced by the second 
moment of the distribution function with the crack orientation vector. A case of special 
interest is a distribution with a rotation symmetry, the uniaxial case. In this case the damage 
tensor can be expressed in terms of a scalar and a unit vector, which is the orientation of the 
rotation symmetry axis. In this case we can define easily a damage vector from the second 
order damage tensor. This damage vector has the orientation of the rotation symmetry axis. 
A representation of the second order damage tensor in the general case without rotation 
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symmetry is also given. It is shown how a vectorial damage parameter can be defined in 
general without rotation symmetry in terms of eigenvectors of the damage tensor. 

2.1. Damage parameter of second order. We define the second order damage tensor as 
the second orientational moment of the distribution function: 



D{x,t) =< I nn >= / I / f{l,n,x,t) nn cfnl^dl 
Jo Js^ 

(2.6) = / [ rm f{l,n,x,t)(fndl 

Jo Js^ 

nn denotes the symmetric traceless part of the dyadic product, and JD is a second order 
symmetric traceless tensor. 

This definition of the second order damage parameter accounts for the crack length 
distribution as well as for the orientation distribution. With an appUed load the length 
distribution evolves in time. This shows up in the evolution of the scalar damage parameter 
as well as in the evolution of the tensor damage parameter. As in our model all cracks rotate 
together with the surrounding material, the time evolution of the orientation distribution is 
a rigid rotation on the unit sphere. 

2.2. Vectorial damage parameter defined from tlie second order tensor. Due to the 

symmetry the second order tensor damage parameter D has a spectral decomposition with 
orthogonal eigenvectors (unit vectors) d, e and /, and eigenvalues 6, e and cj): 



(2.7) D = 5dd + eee + <pff . 
Because D is traceless we have 

(2.8) S + e + (f) = 0. 

Therefore, not all eigenvalues can have the same sign. The following cases concerning the 

signs of the eigenvalues are possible 

(1) One eigenvalue is positive and two eigenvalues are negative, for instance 



(2.9) 



^>0, e<0, $<0 
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In this case we chose the eigenvector (here d) corresponding to the single pos- 
itive eigenvalue as the unit vector defining the orientation of the vector damage 
parameter. 

(2) One eigenvalue is negative and two eigenvalues are positive, for instance 

(2.10) 6>0, e > 0, $ < . 

In this case we chose the eigenvector (here /) corresponding to the single neg- 
ative eigenvalue as the unit vector defining the orientation of the vector damage 
parameter. 

(3) AU eigenvalues are zero: 

(2.11) 6 = 0, e = 0, (t) = 0. 

In this case D = and we have an isotropic orientation distribution. In this case 
no vector damage parameter can be defined. 

(4) One eigenvalue is zero, and the two others have opposite sign, for instance: 

(2.12) 6 = 0, e > 0, = -e < . 

In this very special case we could define a vector damage parameter, having the 
orientation of the wedge-product of the two eigenvectors. 

The length of the damage vector can be defined as the absolute value of the correspond- 
ing eigenvalue. 

The definition of the damage vector in terms of an eigenvector naturally leads to the 
symmetry of an orientation, namely the damage vector and the reversed one cannot be 
distinguished. 

2.3. Special case of uniaxial distribution function. If there exists a rotation symmetry 
axis of the distribution function, two eigenvalues coincide, either the two positive ones, or 
the two negative ones. In both cases the tensor damage parameter is of the form: 

D = {IS{1)) dd = {IS{1)) {dd - ^1 

(2.13) ={lS{l))(^dd-^{dd + ee + ff)^ 

with the unit tensor 1 and a scalar parameter S, denoted as scalar orientational order pa- 
rameter. The unit vector d is the orientation of the rotation synmietry axis. S{1) is a 
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Figure 1 . Different cases of rotation symmetric crack orientation dis- 
tributions, corresponding to different values of the orientational order 
parameter. 

measure of the degree of parallel order of the cracks. It is zero if the orientations are dis- 
tributed isotropically and has the value 1 in the case that all cracks are oriented parallel. 
The orientational order can be different for different crack sizes, therefore S* is a function 
of crack radius I. The average ( ) here is the average over all crack lengths: 

/>oo 

(2.14) {IS{1))^ f{l,x,t)lS{l)fdl. 

Jo 

For the eigenvalues this corresponds to 

(2.15) s^m))(^l-l^^^{lS{l)) 

(2.16) e = -i(/5(0) 

(2.17) ^^_i(Z5(0). 

For positive values of S we have one positive eigenvalue and two negative ones. For 
negative values of S two eigenvalues are positive and one is negative. In both cases the 
definition of the vector damage parameter given in the section above leads to the eigen- 
vector d as the orientation of the damage vector. In case of rotation symmetric orientation 
distributions this is the orientation of the rotation symmetry axis. The case of positive val- 
ues of S corresponds to a distribution, where the crack-normals are more or less parallel to 
the rotation symmetry axis. For negative values of S crack orientations are concentrated in 
a plane perpendicular to the rotation symmetry axis (see figure 1). 

For the damage vector we find in the uniaxial case: 

(2.18) D^^{lS{l))d. 

It depends on the degree of orientational order and on the average crack length. 
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2.4. Case of small deviation of the distribution function from rotation symmetry. If 

the deviation of the orientation distribution from rotation symmetry is small, the two eigen- 
values of equal sign differ only by a small amount s, and the damage tensor is of the form: 
(2.19) 

D = (^-{S{l)l)dd+ [-\{{S{1) - sm)) ee+ {-\{{S{1) - s{l))l)) ff) . 
In this case we can stiU define the damage vector the same way as in the uniaxial case: 

(2.20) 3 ='^{lS{l))d . 

The scalar order parameter 5 is a measure of the degree of order and the biaxiahty param- 
eter s is a measure of the deviation of the orientation distribution from rotation symmetry. 

2.5. Scalar damage parameters. One possible definition of a scalar damage parameter 
is the average crack length: 

POO 

(2.21) D= lf{l,x,t)l^dl . 

Jo 

In the rotation symmetric case another scalar measure of damage is 

POO 

(2.22) Ds = {IS{1)) = / lS{l)f{l,x,t)fdl . 

Jo 

It is the average crack radius, proj ected onto the plane perpendicular to the rotation symme- 
try axis. This is the quantity relevant for describing the progressive damage under external 
load in case of an anisotropic crack distribution. 

3. Examples of constitutive functions for the different damage 

descriptions 

In this section we will show on the example of the free energy density, how the consti- 
tutive equation reduces from a representation in terms of the second order damage tensor 
to one in terms of a vector, or a scalar, respectively. 

We will assume, that constitutive quantities depend on the equilibrium variables strain 
and temperature and in addition on the damage parameter. The temperature dependence 
wiU not be denoted explicitly, as it is a scalar quantity. AH material coefficients may depend 
on temperature. 
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3.1. Free energy as a function of strain and damage in case of a second order damage 
tensor. The most general polynomial form of the energy density up to second order in 
each variable is given by a representation theorem ll3n[32l : 



F{e, D, D) = a2D + y i?^ ^ ^^^^^jj . jj^ 
+ (ai + b^D + C2l?^ + C5tr(£) • £))) tr{e) 
{bs + CiD) tr{e ■ D) + cstr{D D ■ e) 
^+ciD + diD^ + d3tr{D ■ D)^ {tr{e)f 
+ + cqD + daD^ + d4tr{D ■ D)) tr{e ■ e) 
+ (c3tr(e) + d5tr(e • D) + detr{e)D) tr{e ■ D) 
+ (C7 + d^D) tr{e ■ e ■ D) + djtr{e)tr{D D e) 

(3.23) +dgtr{D D e e) 

This form is the simplest and natural extension of linear elasticity considering the dam- 
age. The coefficients still can be arbitrary functions of temperature. 
In the case of rotation symmetry we have: 

(3.24) D = Ds"dd 
and the scalar products can be calculated as: 



D D = Dl dd ■ dd 
(3.25) ^D'sidd- ll) ■ (dd - il) = Dl (Idd + il 



9" 



and 



(3.26) 



tr{D • D) - ^Dl 
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Figure 2. Schematic view of an experiment with uniaxial loading. 



The expression for the free energy simplifies to: 



F{e, Ds, d, D) = a2D + -^D' + -b^Di, 
trie) { ai + b^D + c^D^ - ^b^Dg - ^aDDs + (^05 + ^cs)^^?! 

-{tref ( + ciD + diD^ - ^sDs - ^deDDs 



+i^d3 + ids + ^^7)^1 



11 2 1 

tr{e ■e)[n + ceD + d2D^ - -crDs - -c^DDs + (-^4 + gC^g)^'! 

+e : dd ( b^Ds + c^DDs + ^Cgi^l^ + (e : ddf d^Ds 

+ {tre)e : dd [ [d^ + C'i)Ds + d^DDs + ^d^treDl 



(3.27) 



-d e e d\ cjDs + d^DDs + Jdgi?! 



The free energy density is expressed here in terms of the vector d (the rotation symmetry 
axis) and the scalar damage parameter Ds- 



3.2. Special case of uniaxial strain in the z-direction and symmetry axis of the distri- 
bution in the same direction. This situation occurs (approximately) in a uniaxial tension 
experiment (see figure 2). 

The assumption that in all volume elements the CDF is rotation symmetric with the 
z-direction as symmetry axis is an approximation, valid for small deformations. In case 
of large deformations, the rotation of the volume element cannot be neglected, and the 
cracks rotate with the material element. This leads to a rotation of the symmetry axis of 
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the orientation distribution function in the volume element, which depends on the position 
of the volume element. The (local) symmetry axis of the CDF does not coincide with the 
direction of the applied strain anymore. 

In the geometry with the global rotation symmetry around the z-axis the only interesting 
components of tensors are the z-z-components and traces. For the strain and the damage 
vector we have: 



(3.28) e = e^^e^e^ =: ee^e^ 

(3.29) d = e, 

and the free energy density reduces to a function of three scalar quantities e, D, Ds- D is 
the average crack length (in any direction) and Ds is a measure of the anisotropy of the 
average crack length distribution in the z-direction and orthogonal to it. 

F{e, Ds, D) = a2D + ^D^ + ^b2Dl 
+e (ai + hiD + C2D^ + Q&g + ^CiD^ Ds + Qcs + ^cg^ D^ 

Q + M + (ci + cg) £> + {dl + d2) D'^ 

/2 2 /2 2 

DsUc, + -cr + Dl-d, + -d, 



(3.30) ( -da + -c?4 + -ds + -d^ + -dg ) ) 



2j 2 , 4^ 4 , 4 
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